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We present a high-statistics measurement of differential cross sections and the total cross section 
for the process 77 — » 7r°7r° in the kinematic range 0.6 GeV < W < 4.0 GeV and |cos#*| < 0.8, 
where W and 6* are the energy and pion scattering angle, respectively, in the 77 center-of-mass 
system. Differential cross sections are fitted to obtain information on S, Do, D2, Go and G2 waves. 
The G waves are important above W — f .6 GeV. General behavior of partial waves is studied by 
fitting differential cross sections in a simple parameterization where amplitudes contain resonant 
contributions and smooth background. The D2 wave is dominated by the /2(f270) meson whose 
parameters are consistent with the with the current world averages. The Do wave contains a /2(f 270) 
component, whose fraction is fitted. For the S wave, the /o(980) parameters are found to be 
consistent with the values determined from our recent tv + tt~ data. In addition to the /o(980), the 
S wave prefers to have another resonance-like contribution whose parameters are obtained. 

PACS numbers: 13.60.Le, 13.66.Bc, 14.40.Cs, 14.40.Gx 



I. INTRODUCTION 



Studies of exclusive hadronic final states in two-photon collisions give valuable information on the physics of light and 
heavy-quark resonances, perturbative and nonperturbative quantum chromodynamics (QCD) and hadron-production 
mechanisms. So far. Belle has measured the production cross sections of charged-pion pairs [2, 0, El, charged- and 
neutral- kaon pairs 0, 0, H[ , and proton-antiproton pairs Q . We have also analyzed D- meson-pair production finding 
a new charmonium state identified as the Xc2(2-P) [3]- 

Here we present the cross sections and an analysis of high-statistics neutral-pion pair production in two-photon 
processes. The motivation of this study is essentially the same as that for charged-pion pair production. However, 
the two processes are physically different and independent; we cannot predict one by measuring only the other. 

In the low energy region (W ~ 0.8 GeV, where W is the center-of-mass (cm.) energy of the 7T7t system), it is 
expected that the difference of meson electric charges plays an essential role in the difference between the 7r + 7r~ and 
7r°7r° cross sections. Predictions are not straightforward because of non-perturbative effects. In the intermediate 
energy range (0.8 GeV ~ W ~ 2.4 GeV), formation of meson resonances decaying to 7T7t is the dominant contribution. 
Since ordinary qq mesons conserve isospin in decays to 7T7t, we can restrict the I G J PC quantum numbers of the meson 
produced by two photons to be + (even) ++ , that is, /j= e ven mesons. The ratio of the /-meson's branching fractions, 
&{f ~ * 7r ° 7r °)/'B(/ ~~ * 7r + 7r~) is 1/2 from isospin invariance. However, interference of resonances with the continuum 
component, which cannot be precisely calculated, distorts this ratio even near the resonant peaks. 

A long-standing puzzle in QCD is the existence and structure of low mass (~ 1.5 GeV/c 2 ) scalar mesons Q. They, 
in particular, the /o(600) (a meson), are closely related to the QCD vacuum through spontaneous breakdown of chiral 
symmetry Two-photon-resonant production of a meson gives valuable information such as its two-photon width, 
which is sensitive to its charge structure. Unfortunately, experimental constraints in a B factory experiment do not 
allow measurements much below W ~ 0.8 GeV. 

For higher energies (W ~ 2.4 GeV), we can invoke a quark model. In leading order calculations [l^ which take 
into account spin correlations between quarks, the 7r°7r° cross section is predicted to be much smaller than that of 
7T + 7r~, and the ratio of 7r°7r° to 7r + 7r~ is around 0.03-0.06. However, higher-order or nonperturbative QCD effects 
can modify this ratio. For example, the handbag model, which considers soft hadron exchange, predicts the same 
amplitude for the two processes and hence this ratio is expected to be 0.5 [ll|. Analyses of energy and angular 
distributions of the cross sections are essential to determine properties of the observed resonances and to test the 
validity of QCD models. 
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In this paper, we present measurements of the differential cross sections, do jd\ cos#*|, for the process 77 — > 7r°7r° 
in a wide two-photon cm. energy (W) range from 0.6 to 4.0 GeV, in the cm. angular range, | cos 6 * | < 0.8. The 
95 fb _1 data sample results in several hundred times larger statistics than in previous experiments fl3l]7 The data 
here are concentrated in the low energy region (W < 1.6 GeV). In the higher energy region, many more resonances 
contribute and thus make an analysis based on a single experiment difficult. Furthermore, statistics are still poor in 
this range. The angular coverage up to | cos#*| = 0.8 greatly enhances the capability for separating partial waves. In 
the low energy region, the cross section is dominated by the /2(1270). A clear peak corresponding to the /o(980) is 
found for the first time in two-photon production of 7r°7r°. Furthermore, a "model- independent" partial wave analysis 
(where interference terms of amplitudes are temporarily neglected) reveals another resonance-like structure around 
1.2 GeV in the S wave. This may be due to the contributions of the /2(1270) in the Do wave and/or scalar resonances 
such as the /o(1370). We fit the differential cross sections assuming such contributions and obtain their parameters. 

This paper is organized as follows. In section[TTl the experimental apparatus relevant to this measurement is briefly 
described together with information on the trigger and a description of the event selection. Differential cross sections 
are derived in section ITTT1 In section IIV1 the measured differential cross sections are fitted to obtain the resonance 
parameters of the /o(980) and a scalar meson as well as to extract the fraction of the /2(1270) in the Do wave. Finally 
in section [V] a summary and conclusion are given. 

II. EXPERIMENTAL APPARATUS AND EVENT SELECTION 

We use data that corresponds to an integrated luminosity of 95 fb -1 recorded with the Belle detector at the KEKB 
asymmetric-energy e + e~ collider [T3|. The e + e _ cm. energy of the accelerator was set at 10.58 GeV (83 fb -1 ), 
10.52 GeV (9 fb" 1 ), 10.36 GeV (T(35) runs, 2.9 ib" 1 ) and 10.30 GeV (0.3 ftr 1 ). The differences between the two- 
photon flux (luminosity function) in the measured W regions due to differences in the beam energies are small (at 
most a few percent), and the fraction of integrated luminosity of the runs with lower beam energies is also small; we 
combine the results for different beam energies. The variation of the cross section because of this effect is less than 
0.5%. 

The Belle detector is a magnetic spectrometer covering a large solid angle (polar angles between 17° and 150° and 
the full azimuthal angle). A comprehensive description of the Belle detector is given elsewhere [l5|. We mention 
here only those detector components that are essential to the present measurement. Charged tracks are reconstructed 
from hit information in a central drift chamber (CDC) located in a uniform 1.5 T solenoidal magnetic field. The z 
axis of the detector and the solenoid arc along the positron beam direction, with the positrons moving in the — z 
direction. The CDC measures the longitudinal and transverse-momentum components (along the z axis and in the np 
plane, respectively). Photons are detected and measured in an electromagnetic calorimeter (ECL) located inside the 
solenoid. The ECL is an array of 8736 CsI(Tl) crystals pointing toward the interaction point, which help separating 
photons from 7r°'s for energies up to ~ 4 GeV. 

We require that there be no reconstructed CDC tracks coming from the vicinity of the nominal collision point. 
Photons from decays of two neutral pions are measured in the ECL. Signals from the ECL are used to trigger. The 
ECL trigger requirements are the following: the total ECL energy deposit in the triggerable acceptance region (see 
below) is greater than 1.15 GeV (the "HiE" trigger) or the number of ECL clusters (each crystal having more than 
110 MeV) is four or greater (the "Clst4" trigger). The above energy thresholds are determined from a study of the 
correlations between the two triggers in data. This trigger logic is realized in hardware; no additional software filters 
are applied for events triggered by either of the two ECL triggers. 

The analysis is performed in the "zero-tag" mode, where neither the recoil electron nor the positron is detected. 
We restrict the virtuality of the incident photons to be small by imposing strict transverse-momentum balance with 
respect to the beam axis for the final-state hadronic system. 

The selection conditions for 77 — > ir°ir° signal candidates are the following. All the variables in criteria (l)-(6) are 
measured in the laboratory frame: (1) there is no good track that satisfies dr < 5 cm, dz < 5 cm and p t > 0.1 GeV/c, 
where dr and dz are the radial and axial distances, respectively, of closest approach (as seen in the np plane) to the 
nominal collision point, and the p t is the transverse momentum measured in the laboratory frame with respect to 
the z axis; (2) the events are triggered by either the HiE or Clst4 triggers; (3) there are two or more photons whose 
energies are greater than 100 MeV; (4) there are exactly two r's, each 7r° having a transverse momentum greater 
than 0.15 GeV/c with each of the decay-product photons having an energy greater than 70 MeV; (5) the two photons' 
momenta are recalculated using a 7r°-mass-constraincd fit, and required to have a minimum x 2 value for the fit (there 
was a negligible fraction of events with ambiguous photon combinations); (6) the total energy deposit in the ECL is 
smaller than 5.7 GeV. 

The transverse momentum in the e + e~ cm. frame (|Ep£|) of the two-pion system is then calculated. For further 
analysis, (7) we use events with |Ep£| < 50 MeV/c as the signal candidates. 



3 



In order to reduce uncertainty from the efficiency of the hardware ECL triggers, we set offline selection criteria that 
emulate the hardware trigger conditions as follows: (8) the ECL energy sum within the triggerablc region is greater 
than 1.25 GcV, or all four photons composing the two ir° are contained in the triggerablc acceptance region. Here, 
we define the triggerable acceptance region as the polar-angle (9) range in the laboratory system 17.0° < 9 < 128.7°. 



In this section, we describe the derivation of differential cross sections. First, candidate events are divided into 
bins of W and | cos 9* \ . Backgrounds are then subtracted by fitting the transverse-momentum distribution. Event 
distributions are unfolded to correct for finite energy resolution. Finally, differential cross sections are obtained in 
bins of W and | cos 9* \ . 



We derive the cm. energy W of the two-photon collision from the invariant mass of the two-neutral-pion system. 
We calculate the cosine of the scattering angle of ir° in the 77 cm. frame, | cos#*|. We then approximate the e + e~ 
collision axis in the e + e~ cm. frame as the reference for this polar angle. 

The two-dimensional yield distribution of the selected events is shown as a lego plot in Fig. [T] The W distribution 
with I cos 9*\ < 0.8 is shown in Fig. [5] The total number of events observed is 1.25 x 10 6 . We observe clear peaks for 
the /o(980) near 0.98 GeV and the /2(1270) near 1.25 GeV and find at least two more structures around 1.65 GeV 
and 1.95 GeV. 



FIG. 1: Two-dimensional W and |cos#*| distribution for 7r°7r° candidates in data. The same distribution is viewed from two 
different directions. 



We use the pt -balance distribution, i.e., the distribution in to separate the signal and background components. 

The signal Monte Carlo (MC) shows that the signal component peaks around 10-20 MeV/c in this distribution. In 
the experimental data, however, in addition to the signal component, we find some contributions from pt-unbalanced 
components in the low-VF region. Such pt-unbalanced backgrounds might originate from processes such as 7r 7r°7r°, 
etc. However, the background found in the experimental data is very large only in the low-VF region where the 7r 7r°7r° 
contribution is expected to be much smaller than tt°tt° in two-photon collisions. (Note that a C = — system cannot 
decay to 7r°7r°.) We believe that the backgrounds are dominated by beam-background photons (or neutral pions from 
secondary interactions) or spurious hits in the detector. 

Figures [3] (a) and (b) show the p t -balance distributions in the low W region. With the fit described below, we 
separate the signal components from the background. In the intermediate or higher energy regions, the pt unbalanced 
backgrounds are either less than 1%, buried under the / 2 (1270) peak (Fig. EDJc)), or consistent with zero within 
statistical errors. For the highest energy region 3.6 GeV< W < 4.0 GeV, we subtract a 3% background from the yield 



III. DERIVATION OF DIFFERENTIAL CROSS SECTIONS 



A. Signal Distributions 




B. Background Subtraction 
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FIG. 2: W distribution for candidate events. The angular coverage is | cos#*| < 0.8. The curve is an estimate of backgrounds 
from events with p t imbalance (see text). 

in each bin to account for the background from the pt-unbalanced components and assign a systematic error of the 
same size, although it is not statistically significant even there (Fig. [3jd)). 

A fit to the p t -balance distribution is performed in the region IX^Ptl — 0-2 GeV/c to separate the signal and 
background components for the W region below 1.2 GcV. The fit function is a sum of the signal and background 
components. The signal component is an empirical function reproducing the shape of the signal MC, y = Ax/(x 2A + 
B + Cx), where x = | Y^Pt Ij A B and C are the fitting parameters, and y is the distribution. This function has a 
peak at x ~ (yj)~ and vanishes at x = and as x — * oo. The shape of the background is taken as a linear function 
y = ax for x < 0.05 GeV/c, which is smoothly connected to a quadratic function above x > 0.05 GeV/c. 

The background yields obtained from the fits are fitted to a smooth two-dimensional function of (W 1 |cos0*|), 
in order not to introduce statistical fluctuations. The backgrounds are then subtracted from the experimental yield 
distribution. The background yields integrated over angle are shown in Fig. [2] for W < 1.2 GcV. Above 1.2 GeV, we 
do not find any statistically significant background contributions from the fit. The correction and systematic errors 
from background are summarized in Sect. III.F. We omit the data points in the small-angle (| cos6**| > 0.6) region 
with W < 0.72 GeV, because there the background dominates the yield. 

C. Unfolding the W Distributions 

We estimate the invariant-mass resolutions from studies of the signal MC and data. We find that the MC events 
have a relative invariant-mass resolution of 1.4%, which is almost constant over the entire W region covered by the 
present measurement. The ir° momentum resolution is known to be about 15% worse in the experimental data than 
in MC from a study of the p t -balance distributions. Moreover, the distribution in the MC is asymmetric; it has a 
longer tail on the lower mass side. An asymmetric Gaussian function with standard deviations of 1.9%W / and 1.3%VF 
on the lower and higher sides of the peak, respectively, is used and approximates the smearing reasonably well. This 
invariant-mass resolution is comparable to or larger than the W bin width (20 McV) used in Figs. Q] and [5] We unfold 
the invariant-mass distribution in each | cos 9* | bin separately, to correct for migrations of signal yields to different 
W bins in obtaining the true W distribution, based on the asymmetric Gaussian smearing described above. The 
migration in the | cos 6* | direction is expected to be small and is neglected. 

The unfolding uses the singular value decomposition algorithm [l6[ at the yield level and is applied so as to obtain 
the corrected W distribution in the 0.9 - 2.4 GeV region, using data in the observed W range between 0.72 and 
3.0 GeV. For lower energies, W < 0.9 GeV, migration is expected to be small because of the better invariant-mass 
resolution compared with the bin width. For higher energies, W > 2.4 GeV, where the statistics are relatively low 
and the unfolding would enlarge the errors, we rebin the data with a bin width of 100 McV, instead of unfolding. 
Distributions before and after the unfolding for a typical angular bin (| cos 9* \ — 0.225) are shown in Fig.HJ 

We calibrate the experimental energy scale and invariant-mass distribution using the 77 invariant mass from 
experimental samples of 77' — > 77 in two-photon processes. The peak position is consistent with the nominal mass of 
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FIG. 3: Distribution of imbalance in | 5^P*| f° r candidate events, (a) In the bin centered at W = 0.90 GeV and | cos6>*| = 0.05 
(the bin width is 0.04 GeV and 0.1 in the W and |cos#*| directions, respectively, in (a)-(c)), the experimental distribution 
(dots with error bars) is fitted with the sum of signal and background components (curves). The gray region shows the 
estimated background contamination in the signal region, (b) The same as (a) for the bin centered at W — 0.66 GeV. (c) In 
the W — 1.18 GeV, |cos#*| = 0.65 bin the experimental distribution is compared with the signal MC (histogram), (d) The 
same as (c) for 3.6 GeV < W < 4.0 GeV and | cos6>*| < 0.4. 



rf with an accuracy better than 0.2%. The mass resolution estimated from the peak width is also consistent with the 
smearing function that is used for the unfolding of the 7r°7r° invariant-mass distributions. 



D. Determination of Efficiency 



We determine the efficiency for the signal using the detector and trigger simulators and applying the selection 
criteria to signal MC events. The e + e~ — > e + e~7r°7r° signal MC events are generated using the TREPS code [Tt} 
at 58 fixed W points between 0.5 and 4.5 GeV, and isotropic in | cos#*|. The angular distribution at the generator 
level does not play a role in the efficiency determination, because we calculate the efficiencies separately in each 0.05 
wide | cos &*\ bin. The number of events generated is 4 x 10 5 at each W point. To minimize statistical fluctuations in 
the MC calculation, we fit the efficiency to a two-dimensional empirical function in (W, | cos#*|). The efficiency thus 
determined is depicted in Fig. [5] 

We find that the trigger efficiency, which is defined as the triggerable fraction of events that can pass through the 
event selection criteria, is almost flat and close to 100% for the region W ~ 1.3 GeV as shown in Fig. [5] It decreases 
to typically 30% at the lowest energies, W < 0.7 GeV. Meanwhile, the overall efficiency (the product of the efficiencies 
for the trigger and the acceptance) is about 11% at maximum and decreases (down to around 1%) at lower W or 
smaller cm. angles (larger |cos#*|). 

The overall efficiency calculated using the signal MC events is corrected for a systematic difference found between 
the peak widths in the pt-balance distributions of the experimental data and the MC, which could affect the efficiency 
through the | Y^Pt I cu ^- This originates from a difference in the momentum resolution for 7r°'s between data and MC 
events. We find that the data peak position is 10% to 20% higher than the MC expectation, depending on W and 
| cos#*|. The efficiency correction factor ranges from 0.90 to 0.95. 
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FIG. 4: Invariant-mass distributions before (orange colored triangles) and after (dark-blue diamonds) the unfolding, at | cos 9* | = 
0.225. The inset shows the same distribution on a semilogarithmic scale. The bin width changes from 0.02 GeV to 0.04 GeV 
above W = 1.6 GeV. 



E. Cross-Section Calculation 



The differential cross section for each (W, | cos£>*|) point is derived from the following formula: 

da AY-AB 
d\ cos 9* | ~ AWA\ cos 6»*|J CdtL^W)^ ' ^ ' 

where AY and AB are the signal yield and the estimated p t -unbalanced background in each bin, AW and A| cos6>*| 
are the bin widths, J Cdt and L^IW) are the integrated luminosity and two- photon luminosity function calculated 
by TREPS, respectively, and 77 is the net efficiency. The luminosity function transforms the cross sections for the 
e + e~ incident beam to that of the 77 incident using a relation: 

The W-bin width is 0.02 GeV up to 1.6 GeV and then is modified to 0.04 GeV and 0.1 GeV for W between 1.6 GeV 
and 2.4 GeV and W above 2.4 GeV, respectively. The width of the | cos6>*| bins are fixed to A| cos#*| = 0.05. 

The differential cross sections obtained are shown at several W points in Fig. [71 They show quite different behaviors. 
It should be noted that the cross-section results after the unfolding are no longer independent of each other in 
neighboring bins, in both central values and sizes of errors. 

Figures Ufa) and (b) show the W dependence of the cross section integrated over | cos#*| < 0.8 and | cos#*| < 0.6, 
respectively. They are obtained by adding da/d\ cos#*| • A| cos#*| over the corresponding angular bins. The total 
cross section (| cos#*| < 0.8) is dominated by the /2(1270). A clear peak due to the /o(980) is also visible as in the 
7r + 7r~ cross section from Belle [3, 0- Additional structures are visible around W = 1.65 GeV and W = 1.95 GeV. 
The data points for 0.9 GeV < W < 2.4 GeV are the unfolded results where the bin widths AW are 0.02 GeV 
(0.04 GeV) in W above (below) 1.6 GeV. For the data points above 2.4 GeV, we average five data points each with 
a bin width of 0.02 GeV and obtain results for every 0.1 GeV bin in W. We have removed the bins in the range 
3.3 GeV < W < 3.6 GeV, because we cannot separate the contributions from XcO, Xc2 and the continuum in a 
model-independent way, due to the finite mass resolution and insufficient statistics in this W range. 



F. Systematic Errors in the Cross Sections 

We summarize the evaluation of the systematic errors for a(\ cos6>*| < 0.8) (a(\ cos 8*\ < 0.6) for W < 0.72 GeV) 
at each W point. These come from the following sources: trigger efficiency, reconstruction efficiency p t -balance cut, 
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FIG. 5: Overall efficiency in the W — | cos#*| plane. 



background subtraction, beam-background effect, other efficiency errors including smoothing procedure, unfolding 
procedure and luminosity function. 

Trigger efficiency: the systematic error from the Clst4 trigger is assigned as 2/3 of the difference of the efficiencies 
with different threshold assumptions for the ECL cluster - 110 MeV and 100 MeV - set in the trigger simulator for 
the energy region W < 2.5 GeV. We include a separate 4% uncertainty in the HiE trigger efficiency for the whole W 
region. The systematic errors from the two triggers are added in quadrature. This systematic error becomes large in 
the low W region, 20%-30% for W < 0.8 GeV. 

Reconstruction efficiency: the uncertainty in the ir° reconstruction efficiency is estimated from a comparison of 
£)°-meson decays to K~ir + and K~ir + ir°. An error of 6% for two pions is assigned. 

The pt-balance cut is 3% - 5%, which is one half of the correction discussed above. 

Background subtraction: 20% of the size of the subtracted component is assigned as the error from this source. In 
the W region where the background subtraction is not applied (W > 1.2 GeV), we neglect the error for 1.2 GeV< 
W < 1.5 GeV, and assign 3% for W > 1.5 GeV, which is an upper limit on the background contamination expected 
from the p t -unbalanccd distributions. For W > 3.6 GeV where we have applied a 3% correction for the background 
subtraction, we also assign a systematic error of the same magnitude. 

Beam-background effect for event selection: we assign a 2% - 4% error depending on W for uncertainties of the 
inefficiency in selection due to the effect of beam-background photons. 

Other efficiency errors: an extra error of 4% is assigned for uncertainties in the efficiency determinations based on 
the MC including the smoothing procedure. 

Unfolding procedure: we adopt the change in the unfolded yield in each bin when we modify the effective-rank 
parameter (kset parameter) applied in the unfolding procedure [l6| in a reasonable range as a systematic error. 

Luminosity function: the uncertainty is estimated to be 4% for W < 3.0 GeV and 5% for W > 3.0 GeV. 

The total systematic error is 10% in the wide energy region, 1.0 GeV< W < 3.4 GeV. The error, which is dominated 
by the background subtraction, becomes much larger for lower W, 15% at W = 0.85 GeV, 30% at W = 0.70 GeV 
and 55% at W = 0.61 GeV. For higher W, the systematic error is rather stable, remaining at the 11% level for 
3.4 GeV< W < 4.0 GeV. 
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FIG. 6: Trigger efficiencies plotted as a function of W for two | cos#*| bins. 



G. Comparison of Cross Sections with the Previous Experiment 

The total cross section (|cos#*| < 0.8) is compared with the previous measurement by Crystal Ball at DORIS 
II [HI (Fig. [5]). The agreement is fairly good. The error bars shown are statistical only, and the systematic errors (7% 
for W > 0.8 GeV and 11% for W < 0.8 GeV for the Crystal Ball results) should also be considered in the comparison. 
The present measurement has several hundred times more statistics than the Crystal Ball measurement. 

IV. FITTING DIFFERENTIAL CROSS SECTIONS 

In this section, differential cross sections at each W bin are fitted to obtain information on the partial waves. A 
simple parameterization is then fit to differential cross sections to obtain resonance parameters of the /o(980) and 
another scalar meson denoted here as fo(Y) as well as to extract the /2(1270) fraction in the Do wave. 

A. Formalism 

In this channel, only partial waves of even angular momenta contribute. Furthermore, in the energy region W ~ 
3 GeV, J > 4 partial waves (next is J = 6) may be neglected so that only S, D and G waves are to be considered. 
The differential cross section can then be expressed as 



da 
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D 2 Yi + G 2 Yl 
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FIG. 7: Differential cross sections for the eight selected W points indicated in GeV in each figure. The results for W < 2.4 GeV 
are after the unfolding. The curves are fit results described in sec. IIVI the solid (dashed) lines are the result of the SD (SDG) 
fit (see text), while dotted lines indicate the size of systematic errors in the cross sections. 



|cos 8*| < O.i 



(a) 



|cos e*| < 0.6 



(b) 
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FIG. 8: The total cross section for (a) | cos#*| < 0.8 and (b) | cos#*| < 0.6 after the unfolding and rebinning. See text for the 
bins in the range 3.3 GeV < W < 3.6 GeV. 
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FIG. 9: Cross section integrated over the angular region | cos 9* | < 0.8 compared with the previous measurement from the 
Crystal Ball [l3j]. Dashed lines show the total systematic errors for the Belle measurement. Systematic errors of similar size 
(not shown) are present in the Crystal Ball measurement. 

where Do and Go (Aa and G2) denote the helicity (2) components of the D and G waves, respectively, and YJ 1 are 
the spherical harmonics: 



1 



Y° 
J 2 



-^-(3cos 2 r-l), 

l07T 



\Y4 



Y° 



15 
32^ 



■ sin 



— J- (35 cos 4 6* - 30 cos 2 6* +3) 

16 V 7T 



(7 cos 2 6* - 1) sin 2 9* 



8 V 2tt 



(3) 



Since the I^J^Is are not independent, partial waves cannot be separated using measurements of differential cross 
sections alone. 

To overcome this problem, we begin by rewriting Eq. ([2]) as: 



da 



( 77 - tt'v-) = s 2 \Y r + m i^t + m m + ^ mT + g\ 



(4) 



Awd\ cos 6 

The amplitudes S 2 , etc. correspond to the cases where interference terms are neglected. When interference terms are 
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included, they can be written as (see Appendix IX]) : 

S 2 = \S\ 2 + V5di(S*D Q ) - 45R(S*G ) + -^=9t(D^G ) + ^^^{D* 2 G 2 ) , 

V 5 5 

Dl = \D Q \ 2 + ^(S*D ) + 2^(S*Go) + ^(D* Go)-^-^(D* 2 G 2 ) , 

D 2 = \D 2 \ 2 -^&{S*Do) + 2&{S*G Q ) + V5M(DZG )-^$t(DZG 2 ) , 

G 2 = \G Q \ 2 + -j=n(D*G ) + ^n(D* 2 G 2 ), 

G\ = \G 2 \ 2 - V5ft(D*G ) - -j=n(D* 2 G 2 ) . (5) 

Since squares of spherical harmonics are independent of each other, we can fit differential cross sections at each W to 
obtain S 2 , D 2 , £>f, G 2 , and G\. The fit up to J = 4 is called the "SDG fit". At low energy, we expect that J = 4 
waves are unimportant. Therefore we also perform a separate fit setting G 2 = G\ = 0, which is called the "SD fit". 

The unfolded differential cross sections are fitted, where only statistical errors are taken into account in the fit. 
Although they are not independent at each W because of the unfolding procedure, we treat them as independent in the 
fit. The effect of correlations between bins is taken into account in systematic errors as described below. Differential 
cross sections for | cos#*| < 0.8 are available for 0.73 GeV < W < 4.0 GcV. 

Examples of the fit quality are shown in Fig. [7] In the energy region near W = 2 GeV, a need for J = 4 waves is 
evident. There, G\ deviates from zero as can be seen in Fig. [TU] Since the behaviors of \Y®\ 2 and \Y 2 \ 2 are rather 

similar for | cos#*| ~ 0.7, we also fit with |F 4 °| 2 ±|F 4 2 | 2 . The bump in Gq+G 2 . may indicate the presence of the / 4 (2050). 
However, in the high energy region W > 1.6 GcV, there are many more resonances contributing and thus the model 
uncertainty becomes much larger. Therefore in this paper we focus on the energy region 0.8 GcV < W < 1.6 GcV, 
where G waves can be neglected. The S* 2 , Dq and D 2 spectra are shown in Fig. ITTI 

B. Fitting Partial Wave Amplitudes 

In this section, we derive some information on the relevant resonances by fitting differential cross sections by 
assuming certain functional forms for the partial wave amplitudes to understand the general behavior of partial wave 
amplitudes and to check the consistency with the 7r + 7r~ data [H, 0. Note that we do not fit the obtained S 2 , etc. but 
fit the differential cross sections directly; once the functional forms of amplitudes are assumed, we can use Eq. J2]) to 
fit differential cross sections. 

Here, we concentrate on the energy region W < 1.6 GeV, where Go and G2 can be neglected. A full amplitude 
analysis will be performed in the near future using all data available in addition to the 7r 7r° data. Thus we employ 
a simple parameterization in this study. The D 2 spectrum is dominated by the /2(1270) resonance. The /2(1270) 
could contribute to the Do wave. The S 2 distribution has peaks apparently due to the /o(980) resonance and another 
resonance-like structure around 1.2 GcV, which motivates us to introduce a scalar meson denoted here as fo(Y) 
together with a contribution from the /2(1270) in the Do wave. The fo(Y) can be either the /o(1370) or /q(1500) or a 
mixture of both [l8[. Note that the mass and width of the /o(1370) are known with a large uncertainty [llj]. Neither 
of these states has been observed previously in two-photon production [l9[. The goal of this analysis is to obtain 
parameters of the /o(980) and fo(Y), to check the consistency of the /2(1270) parameterization where its two-photon 
branching fraction is floated and to measure the helicity 0-to-helicity 2 ratio of /2(1270) production. 

1. Parameterization of Amplitudes 

Based on the above observations, the S, Do and D2 waves are parameterized as follows: 

S = A M980) e^>+A MY) e^+B s , 



Do = J-^A Ml270) e^+B D0 , 
u 1 + r 02 



D 2 = \/ Y^ A f2(i270)e^ d2 + A mi525) e^' + B D2 , (6) 
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where A^ ( 980 ), Af ( Y ), -<4f 2 (i27o) an( i ^./ 2 (i525) are the amplitudes of the /o(980), another scalar resonance denoted as 
fo(Y), the /2(1270) and the /2(1525), respectively; Bs, Boo and Bp 2 are "background" amplitudes for S, D and D 2 
waves; and <f> S Q, <p s \, 4>d0i 4>d2 an d 4>d2' are the phases of resonances relative to background amplitudes. The parameter 
ro 2 represents the fraction of the / 2 (1270) component in the Do wave. Here, as a default, we assume the presence of 
both the fo(Y) in the S wave and the / 2 (1270) in the Do wave. We also study the cases where either ro2 = or there 
is no f (Y). 

We assume background amplitudes to be quadratic in W for all the waves: 

B s = a sr W 2 + b sr W + c sr + i{a sl W 2 + b sl W + c sl ), 
B D0 = a W 2 + b Q W + c , 

B D2 = a 2 W 2 + b 2 W + c 2 . (7) 

The D and D 2 background amplitudes arc taken to be real by definition. 

We use the parameterization of the /o(980), /2(1270) and /2(1525) given in Ref. [l| and 0]. We note that 
B(fj — > 7r°7T )/$(/,/ — > 7r + 7r _ ) = 1/2 (because the fj mesons are isoscalars). For completeness, we reproduce 
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FIG. 11: The S 2 (top), Dq (middle) and D 2 (bottom) spectra. Curves are result of the fit with the parameterization discussed 
in text. The error bars are diagonal statistical errors only. 



the parameterization of the /o(980) and /2(1270). For the /o(980) meson, we take 



Vo(980) 



i 



w 



16tt 



D 



(8) 



where (3x = y 1 — 4 ~^~2~ is the velocity of the particle X with mass Mx in the two-body final states, gj xx is related 
to the partial width of the /o(980) meson via Txx{fo) = \^m~, X ■ ^ ne f ac tor Df is given as follows [IT 



D fo (w) = m% — w 2 + m£ (m /o ) - n£> (w) + 5?n£ (M fo ) ~ n£ (w) , 



where for X = ir or K, 



11% (W) 



Px9 2 fo xx 



16-k 



'., 1 , l-Px 
i-\ In 

7T l + px 



(9) 
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The factor (3 k is real in the region W > 2Mk and becomes imaginary for W < 2Mk- The mass difference between 
K± and K° (K°) is included by using (3k = \ (/?if± + Pk° ) ■ The parameters assumed and determined in Ref . [l| are 
summarized in Table [U 



TABLE I: Parameters of the ft (980) assumed and fitted in Ref. Q. 

Parameter Value Unit Reference 

Mais' 985.6 ±\l (stat) ±\\ (syst) MeV/c 2 [1] 
5/o(980)„ 1-33 (stat) tUt (sy^) GeV [1] 
r^,- 34.2 t\H (stat) t*; 8 5 (syst) MeV [1] 
5/o(980)xk/3/ (980)^ 4.21 ±0.25 (stat) ± 0.21 (syst) [21J 
r%7 205 ±|g (stat) (syst) eV [1] 



Next, we give the parameterizations of the /2(1270) and /2(1525) mesons. The relativistic Breit-Wigner resonance 
amplitude Ar(W) for a spin- J resonance R of mass rriR is given by 



- V W X m\-W^-im R YUW) ' (10) 

Hereafter we consider the case J = 2 (the /2(1270) and /2(1525) mesons). The energy-dependent total width Y tQt (W) 
is given by 

T tot (W) = J2 T xx(W) , (11) 

x 

where X is n, K, 7, etc. The partial width T x x (W) IS parameterized as [Hf 

T x M-WR~XK(*m)' *^'\ . ,12, 

\1x(m R )J D 2 {qx{m R )r R ) 



where T R is the total width at the resonance mass, qx(W 2 ) = \JW 2 /4 — m x , D 2 (x) = 1/(9 + 3x 2 ±.t 4 ), and r R is an 
effective interaction radius that varies from 1 (GeV/c) -1 to 7 (GeV/c) -1 in different hadronic reactions [2j|. For the 
47r and the other decay modes, T^iW) = TrB(R — ► 47r)^- is used instead of Eq. (fT2]l . In Ref. @, all parameters of 

the /2(1270) are fixed at the PDG values [l|j], except for tr, as summarized in Table [TT1 

In the fit below, we float the branching fraction of the /2(1270) into two photons, because its value determined in 
the past experiments is based on various assumptions. 



TABLE II: Parameters of the ft (1270) and ft (1525) assumed 

Parameter / 2 (1270) ft (1525) Unit Reference 

Mass 1275.4 ± 1.1 1525 ± 5 MeV/c 2 [18] 

r to t 185.2±1;| 73± 6 5 MeV [18] 

B{ft^nn) 84:.8±H 0.82 ±0.15 % [18] 

B{ft->KK) 4.6 ±0.4 88.8 ±3.1 % [18] 

B{ft -> m) 10.3 ± 3.1 % [18] 

B(ft^jj) 14.1 ±1.3 1.11 ±0.14 xlO -6 [18] 

r fl 3.62 ±0.03 (GeV/c)" 1 [2] 



Finally, the parameterization of the fo{Y) meson is taken to be: 



/ 87rm /o(y) yr /o(y) r 77 (/ o(y))g(/ (y) -> 

where T 77 (/o(F)) is the two-photon width of the fo(Y) meson. 



AfoiY) ~ m l iY) -w*-im MY) r MY) ' ^ 
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FIG. 12: Total cross section (|cos#*| < 0.8) and results of parameterization. Contributions of S 2 , £>q and t)\ are also shown. 



TABLE III: Fitted parameters 



Parameter 


Nominal 


ro2 = 


No f (Y) 


Unit 


Mass(/ (980)) 
F 77 (/ (980)) 


982.2 ± 1.0 

285.5±i?:? 


980.2 ±1.0 
297.0^37 


983.7+iJ5 
370.5± 2 £ 2 . 


MeV/c 2 
eV 


3/o(980)7T7r 


1.82 ±0.03 


1.79 ± 0.03 


1.89 ±0.03 


GeV 


Mass(/ (Y)) 


1469.7 ± 4.7 


1466.8 ± 0.6 




MeV/c 2 


r(/o00) 

r 77 S(/o(F)^7r°7r°) 


89-7±S"i 
11 2+ 50 

1J -- Z — 4.0 


422.4±i|;| 
6780.2t^6.s 


(fixed) 


MeV 
eV 


ro2 

B(/ 2 (1270)^77) 


o.uy_ 29 

1.57 ±0.01 


(fixed) 

1 62+ 002 


5 04 +u -^ 

, r 9 +0.13 


% 
xlO" 5 


X" (ndf) 


1010.1 (615) 1206.1 (617) 1253.3 (619) 





2. Fitted Parameters 

We fit differential cross sections with the parameterized amplitudes in the cm. energy range 0.8 GeV< W < 
1.6 GeV. The parameter g\ K i 'g^ for the /o(980) is fixed to zero because this is preferred by the fit with a large error. 
There are 25 parameters to be fitted. The coefficients ao and 02 in the continuum paramctcrizations arc chosen to 
be positive to fix sign ambiguities. About 3000 sets of randomly generated initial parameters are prepared and fitted 
using MINUIT [24| to search for the true minimum and to find any multiple solutions. Once solutions are found, 
several tens of MINUIT iterations are needed for convergence; with many parameters (25 here and 21 — 26 for later 
analyses); the approach to the minimum is rather slow. A unique solution is found with 'x 2 /ndf = 1010.1/615 = 1.64 
for the nominal fit, where ndf is the number of degrees of freedom. The fitted parameters are listed in Table Hill The 
errors quoted are statistical only. The fit quality is adequate, x 1 /ndf = 1.64, and represents the trend of the squared 
amplitudes as shown in Fig. 1111 The error bars in the figure are diagonal statistical errors only. The quantity D\ 
is well reproduced except below 1.1 GeV. The effect of the /2(1525) is rather small and not visible in the figures. 
Additional assumptions or a more complicated model are needed to better reproduce the structures visible in Dq for 
the range 1.1 GeV ~ W ~ 1.4 GeV (the / 2 (1270) region). 
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The total cross section (| cos#*| < 0.8) can be obtained by integrating Eq. (|4|) as 

cr tot = 0.8S 2 + 0.457281>o + 0.988288Z>2 , (14) 

where the numerical factors come from the integration of spherical harmonics for | cos6>*| < 0.8. The measured total 
cross section is compared with the prediction obtained from the fitted amplitudes as shown in Fig. [12] They are 
reproduced reasonably well. 

So far we assumed a need for both the fo(Y) in the S wave and the /2(1270) in the Do wave (ro2 ^ 0). We also 
study cases where either one of them is absent. One thousand sets of randomly generated initial values are prepared 
and fitted to find the true minimum. Unique solutions are found in each case. The values obtained are also listed in 
Table IHT1 The \ 2 /ndf clearly favors the presence of both components. 

Comparing the case of no fo(Y) and no /2(1270) in the Do wave, the fit prefers the latter: \ 2 /ndf = 1206.1/617 
compared to 1253.3/619. Thus we conclude that the possibility of only the /2(1270) in the Do wave is disfavored 
compared to the case of only the fo(Y) in the S wave by more than 6 standard deviations, which is calculated from 
the difference of the \ 2 taking the difference of ndfs into account. 

3. Study of Systematic Errors 

The following sources of systematic errors on the parameters are considered: dependence on the fitted region, nor- 
malization errors in the differential cross sections, assumptions on the background amplitudes, uncertainties from the 
unfolding procedure, uncertainties in the parameters assumed for the /o(980) and uncertainties in the measurements 
of the / 2 (1270) and $(1525). 

In each study, a fit is made allowing all the parameters to float and the differences of the fitted parameters from 
the nominal values are quoted as systematic errors. One thousand sets of randomly generated input parameters are 
again prepared for each study and fitted to search for the true minimum and for possible multiple solutions. Unique 
solutions are found in each case. 

Two fitting regions are tried: a narrow one (0.84 GeV < W < 1.56 GeV) and a wide one (0.76 GeV < W < 
1.64 GeV). Normalization error studies are divided into those from uncertainties of the overall normalization and those 
from smearing of the spectra in either | cos#*| or W. For overall normalization errors, fits are made with two sets of 
values of differential cross sections obtained by multiplying by (l±a e (w\ CO s0*|))> where a e is the relative efficiency error. 
For smearing studies, ±4% errors are assigned and differential cross sections are smeared by (1 ± . 1 1 cos0*| =F 0.04) 
and (1±0.W=F0.12). 

For studies of background (BG) amplitudes, one of the waves is changed to either a first- or a third-order polynomial. 
In estimating uncertainties from the unfolding procedure, the kset parameter [l6[ is varied by ±4 from its nominal 
value. When this parameter is increased, the constraints between adjacent bins become weaker and oscillating solutions 
tend to appear as the statistical errors increase. When kset decreases, the opposite trend occurs. 

For the parameterization of the /o(980), the fit prefers zero for the ratio g 2 KK i 'g 2 ^ as stated above and hence we set 
the ratio to be zero in the nominal fit. A systematic error due to this assumption is studied by setting the ratio to be 
4.45, a value 1 standard deviation larger according to the BES measurement [2l|. Finally, the /2(1270) and $(1525) 
parameters are varied by their errors. 

The resulting systematic errors are summarized in Table HVl It is noted that the mass of the fo(Y) jumps from 
~ 1.5 GeV/c 2 to ~ 1.2 GeV/c 2 , which is "compensated" by an increased value of r 2 in a systematic error study 
on unfolding. There, the change of cross-section values from the nominal ones is rather small (at most 15% at 
| cos 6*\ ~ 0.8) with some systematic dependence on cos 9* in some W bins, showing the sensitive nature of this kind 
of analysis. Total systematic errors are calculated by adding individual errors in quadrature. 

4- Summary of Fit Results 

Once the amplitudes are parameterized, differential cross sections can be fitted to obtain the parameters as described 
above. The results are much more powerful than a simple fit to the total cross section. This is because there are so 
many points available that provide rich information. Although the fit quality is not very good as can be seen from 
X 2 /ndf = 1.64, the fit is stable despite the fact that the approach to the minimum is slow requiring tens of MINUET 
iterations. 

In Tables IVl and |VT1 the results obtained for /o(980) and fo(Y) ar e summarized and compared with the PDG (l8j 
and with previous measurements [ U_j25l |. The /o(980) parameters obtained here are consistent with those obtained 
from Belle's measurement of tt + tt~ [l|. The fo{Y) fitted mass is close to the /o(1500) mass, but is also consistent with 
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TABLE IV: Systematic uncertainties for /o and parameters 





/o(980) 


MY) 






Source 




r 


.y7T7r 


Mass 


J- tot 


r B n „ 


n>2 


B(/ 2 -> 77) 




(MeV/c 2 ) 


feVI 




(MeV/c 


2 ) (MeV) 
) \+ vj -° v / 


feV) 


(%) 


(xl0~ 5 ) 


VK-range 


+ 0.0 


+ 11S.8 


+ 0.10 


+30.0 


+ 1.4 


+ Y.5 


+U.51 


+0.57 


-2.8 


-0.0 


-0.02 


-49.7 


-13.0 


-5.0 


—0.00 


-0.00 


Normalization 


+0.6 
—0.4 


+0.0 

— 28.5 


+0.00 
— 0.01 


+8.8 
—0.0 


+ 15.0 
—3.2 


+10.5 
— 1.5 


+0.29 
— 1 61 


+0.44 
— 14 


Bias:| cos 9* \ 


+0.0 


+6.9 


+0.00 


+ 16.0 


+ 13.1 


+8.5 


+0.00 


+0.38 


— 1.0 


— 28.9 


— 0.02 


— 0.0 


—0.0 


— 0.0 


-1.09 


-0.00 


B\w,:W 


+0.1 


+ 10.8 


+0.00 


+25.8 


+ 1.0 


+8.2 


+0.01 


+ 0.58 


— 2.2 


-0.0 


-0.01 


-2.5 


— 2.4 


-0.3 


-0.15 


-0.01 


BG-M(S) 


+0.0 

-1.3 


+ 127.1 
-9.3 


±0.02 


+6.8 
-19.5 


+0.3 
-3.1 


+2.2 

-3.8 


+ 1.14 
— 1 66 


+0.07 

— 01 


BG:S(S) 


+0.0 
-2.5 


+35.3 

-0.0 


+0.00 
-0.01 


+ 19.8 

-0.0 


+0.0 
-7.6 


+0.8 
-0.6 


+2.52 
—0 00 


+0.02 
— 0.01 


BG:D 


+0.1 


+30.6 


+0.01 


+27.1 


+ 1.3 


+8.4 


+0.26 


+0.59 


— 2.6 


— 5.3 


— 0.00 


— 2.2 


— 1.7 


— 0.5 


—0 62 


— 00 


BG:D 2 


+0.0 


+30.7 


+0.01 


+42.6 


+6.3 


+ 10.0 


+0.00 


+0.51 


— 2.9 


— 0.0 


— 0.00 


— 0.0 


— 10.6 


— 0.0 


-0.38 


-0.00 


Unfolding 


+0.0 


+ 2.4 


+0.22 


+0.0 


+39.6 


+602.9 


+ 15.63 


+0.33 


-5.0 


-48.5 


-0.16 


-249.5 


-0.0 


-0.0 


-0.37 


-0.00 


9KK 


+8.0 


+97.4 


+0.00 


+ 10.5 


+ 12.7 


+ 11.1 


+0.00 


+ 0.46 


-0.0 


-0.0 


-0.02 


-0.0 


-0.0 


-0.0 
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/2:mass 


+0.1 


+ 26.5 


+0.02 


+5.1 


+ 12.3 
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+0.36 


+0.00 


-0.7 


-26.6 


-0.03 


-3.8 


-10.7 


-2.9 


—0 21 


-0 01 


ji :width 


+0.6 
-1.2 


+3.9 
-2.8 


+0.00 
-0.01 


+11.7 
-10.3 


+2.6 
-0.5 


+ 1.5 

-0.0 


+0.13 
—0 11 


+0.01 
— 02 


H ■ B(7T7T) 


+0.0 
—0.1 


+0.7 
— 0.0 


±0.00 


+0.4 
— 0.0 


+ 1.7 
— 0.0 


+0.7 
— 0.0 


+0.00 
—0 07 


+0.02 
— 0.00 


h ■ B{KK) 


+0.0 
—0.3 


+4.0 
— 8.7 


+0.00 
— 0.01 


+1.8 
— 0.0 


+2.5 
—0.0 


+ 1.4 

-0.1 


+0.02 
-0.00 


±0.00 


h ■ tr 


±0.1 


+26.5 
-3.3 


+0.02 
-0.01 


+5.1 
-0.0 


+ 12.3 
-0.0 


+ 5.7 
-0.0 


+0.36 
-0.08 


±0.00 


/^mass 


+0.0 
-0.1 


+0.9 

-0.5 


+0.00 
-0.00 


+ 1.0 

-0.6 


+0.3 
-0.0 


+0.3 
-0.1 


+0.05 
-0.00 


±0.00 


/2 :width 


±0.0 


±0.5 


±0.00 


+0.5 
-0.3 


+1.1 

-0.4 


+0.2 
-0.0 


±0.00 


±0.00 


fi ■ B(nn) 


±0.0 


+ 1.3 

-0.8 


±0.00 


+2.3 
-1.9 


+0.2 
-0.9 


+0.9 
-0.8 


+0.05 
-0.02 


±0.00 


f2 ■ B{KK) 


±0.0 


+ 0.5 
-0.0 


±0.00 


±0.1 


+0.2 
-0.0 


+0.1 
-0.0 


±0.00 


±0.00 


fi ■■ Birr) 


±0.0 


+0.7 
-0.2 


±0.00 


±1.4 


+0.6 
-0.7 


+0.7 
-0.6 


±0.01 


±0.00 




±0.0 


+0.2 
-0.1 


±0.00 


±0.0 


+0.1 
-0.0 


±0.0 


±0.00 


±0.00 


Total 


+8.1 


+ 210.9 


+0.24 


+V2.1 


+49.9 


+603.4 


+15. «y 


+1.39 


-8.0 


-70.1 


-0.17 


-255.4 


-22.0 


-7.2 


-2.85 


-0.14 



the /o(1370) mass because of the large systematic error in this experiment and the large uncertainty in the /o(1370) 
mass from the PDG [H|. Although the product r 77 fi(7r°7r ) is consistent with zero given the size of the systematic 
error, the possibility that r 77 2?(7r 7r°) = is disfavored according to the fit (by comparing the nominal fit and the fit 
with no f (Y) in Table III). 

The branching fraction of the /2(1270) to two photons is measured to be (1.57±0.01 to'u) x 10 -5 m good agreement 
with the value (1.41±0.13) x 1(T 5 in the PDG [3]. The value of r 02 , the helicity 0-to-helicity 2 ratio of the / 2 (1270), 
is r 2 = (3.7 ± 0.3 1~2 5 9 9 )%- This is the first measurement that does not neglect interference. However, the large 
systematic errors in the fo{Y) parameters and ro2 indicate the subtle nature of this kind of fitting. We find that the 
parameters of the fo(Y) and the /2(1270) in Do shown in Table III are rather strongly correlated. There seem to 
be structures that require helicity=0 components both around 1.2 GeV and 1.4 GeV. The structure near 1.2 GeV 
is more prominent, which is supported by the relatively robust ro2 value. When a different unfolding solution is 
fitted, the mass of the fo{Y) jumps from 1.4 GeV to 1.2 GeV due to small changes in the cross sections and their 
statistical errors. It is difficult to disentangle these two structures just based on a simple model. We simply quote 
large systematic errors for the fo(Y) parameters. 



TABLE V: Fitted parameters of the /o(980) 



Parameter 


Belle(7r ,J 7r u ) Belle(7r+7r") PDG Unit 


Mass 

1 77 


982.2 ± 1.0 +| o 985.6 +L { 2 5 + ^ 980 ± 10 MeV/c 2 
286±17i? 1 205 li i 1 ^ 310 ±*° eV 
66.9 ± 2.2 tlH 5L3 ii?!? ±3 3 8 2 Unknown MeV 
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TABLE VI: Fitted parameters of the f (Y) 



Parameter 


Belle(7r u 7r u ) Crystal Ball / (1370)(PDG) / (1500)(PDG) Unit 


Mass 
Ttot 
r 77 B(7r°7T ) 


1470 ±% t ass 1250 1200 - 1500 1507 ± 5 MeV/c 2 
90 t\ tl° 2 268 ± 70 150 - 200 109 ± 7 MeV 
11 t| tf 3 430 ± 80 Unknown Not seen eV 



V. SUMMARY AND CONCLUSION 

We present the total and differential cross sections for the process 77 — ► 7r°7r° for 0.6 GeV < W < 4.0 GeV with the 
Belle detector at the KEKB asymmetric-energy e + e~ collider. The 95 fb _1 data sample has several hundred times 
higher statistics than the previous measurements. The differential cross sections are measured up to |cos#*| = 0.8, 
which gives high sensitivity to the behavior of amplitudes. A clear peak corresponding to the /o(980) is observed 
besides the dominant /2(1270) and a dip-peak structure around W = 1.6 GeV in the total cross section. A general 
behavior of amplitudes is studied by fitting the differential cross sections in a simple model, which includes the S, Do 
and D2 waves that are parameterized as smooth backgrounds and resonances: the /o(980), /2(1270), /2(1525) and 
another possible scalar resonance. 

We obtain a reasonable fit with the /2(1270) parameters fixed at the world-average values and its two-photon width 
floating. The fit yields its value, which is consistent with the world-average one. The /o(980) parameters fitted are 
consistent with the values determined in the 7r + 7r~ analysis [H, 0]- Note that the latter are obtained just by fitting the 
total cross section, while in this paper we fit differential cross sections. The structure in S* 2 around 1.2 GeV can be 
reproduced by the fraction of the /2(1270) present in the D wave and/or the fo(Y). However, we cannot disentangle 
them more clearly. 
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APPENDIX A: INTERFERING AMPLITUDES 

In order to obtain Eq. ©, we express F 2 °, F 4 °, YJ?Y£ and \Y?\\Y£ |, in terms of |F °| 2 , |F 2 °| 2 , |F 2 2 | 2 , |V 4 °| 2 , and |F 4 2 | 2 . 
We use a power series in cos 2 8*. The functions, \Y$\ 2 , \Y£\ 2 , {Y^] 2 , \Y°\ 2 , and \Y 2 \ 2 can be written as 

/4rr|y °| 2 \ 
4tt\Y 2 °\ 2 
4tt\Y 2 \ 2 
47r|y 4 °| 2 
\47r|y 4 2 | 2 / 



/ 1 
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■x 

81 
I.I 
45 
.52 
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2211'. 
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cos 4 e* 

cos 6 e* 
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Equation (IA1[) is inverted to obtain 



cos" <r 
cos 6 9* 
V cos 8 e* 



( 1 
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15 
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15 


16 


73 
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V 175 1225 1225 












2 

15 








4 
15 








34 


16 


8 


175 


1575 


315 


184 


256 


32 


1225 


99225 


735 



/47r|y °| 2 \ 



0|2 

4ir\Y 2 °\ 2 



(A2) 



47r|Y 4 °| 2 

V47r|y 4 2 



I 2 / 



Equation ((5|) can then be derived. 
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